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ABSTRACT 

The central region of the Milky Way Galaxy provides a unique laboratory 
for a systematic, spatially-resolved population study of evolved massive stars of 
various types in a relatively high metallicity environment. We have conducted a 
multi-wavelength data analysis of 180 such stars or candidates, most of which 
were drawn from a recent large-scale HST/ NICMOS narrow-band Paschen-a 
survey, plus additional 14 Wolf-Rayet stars identified in earlier ground-based 
spectroscopic observations of the same field. The multi-wavelength data include 
broad-band infrared (IR) photometry measurements from HST/NIC2, SIRIUS, 
2MASS, Spitzer/IRAC, and X-ray observations from Chandra. We correct for 
extinctions toward individual stars, improve the Paschen-a line equivalent width 
measurements, quantify the substantial mid-IR dust emission associated with 
WC stars, and find X-ray counterparts. In the process, we identify 10 fore- 
ground sources, some of which may be nearby cataclysmic variables. The WN 
stars in the Arches and Central clusters show correlations between the Paschen-a 
equivalent width and the adjacent continuum emission. However, the WN stars 
in the latter cluster are systematically dimmer than those in the Arches cluster, 
presumably due to the different ages of the two clusters. In the equivalent width- 
magnitude plot, WNL stars, WC stars and OB supergiants roughly fall into three 
distinct regions. We estimate that the dust mass associated with individual WC 
stars in the Quintuplet cluster can reach 1O~^M0, or more than one order of 
magnitude larger than previous estimates. Thus WC stars could be a significant 
source of dust in the galaxies of the early universe. Nearly half of the evolved 
massive stars in the Galactic center are located outside the three known massive 
stellar clusters. The ionization of several compact HII regions can be accounted 
for by their enclosed individual evolved massive stars, which thus likely formed 
in isolation or in small groups. 
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1. Introduction 

Massive stars play an important role in galaxy formation and evolution. They are 
a primary source of ionizing photons, radiation pressure, mechanical energy, and chemical 
enrichment of the interstellar medium (ISM). When a massive star reaches its evolved stage 
(off the main sequence; MS hereafter), such feedback has its greatest impact, profoundly 
shaping the local and possibly even galactic environment. 

Stellar evolution theories, though with significant uncertainties, suggest that there are 
several phases through which a very massive star with an initial mass > 25Mq should evolve 
between leaving the MS and exploding as a supernova. Depending on initial stellar mass 
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Crowther! 120071 ) . WN and WC represent two main classes of Wolf Rayet (WR) stars: 
the spectrum of a WN star is characterized by strong emission lines of nitrogen (NIII — V) 
relative to Hel — II, whereas a WC star exhibits strong carbon lines (CIII — IV). These two 
classes are further divided into early (WN2-5: WNE and WC4-6: WCE) and late (WN7-9: 
WNL and WC7-9: WCL) types. In particular, the WCL and WNL types seem to be favored 
in a high metallicity environment (jHadfield et al.l 120051 ). The emission lines are believed 
to arise from optically-thin stellar winds, which can also produce a flat free-free continuum 



strong enough to flatte n the stellar spectral shape in the mid-IR ( 1 Wright fc Barlowlll975 



Mauerhan et al. 



stars (iCrowther 



2010a[) . In addition, significant dust emission is often found around WC L 



Tuthill et al.ll2006h . 



20031 ) and makes them very red in the mid-IR (e.g 
The spectra of Ofpe/WN9 stars are similar to those of OB supergiants, but show strong 
HeII(4696A) emission lines. Stars in the LBV phase experience sporadic, giant eruptive ac- 
tivities and are much rarer than WR stars. Thus the LBV phase should be very short. Much 
is yet to be learned to quantify the strong feedback from these classes of evolved massive 
stars (EMSs) as well as their evolution and relative populations. 

Existing studies of EMSs are very much limited by small sample sizes. Because of 
their short lifetimes, EMSs are found abundan t only in starburst re gions. However, even in 
nearby nuclear starburst galaxies such as M83 ( iHadfield et al.ll2005[ ). such regions cannot be 
adequately resolved to identify individual EMSs, especially if they are in compact massive 
stellar clusters. Local Group galaxies are mostly in relatively quiescent states; large samples 
of EMSs have to be collected from diverse regions having different star formation histories 
and environments. Such samples are not suitable for a statistical study of EMSs (e.g., 
to determine their evolutionary dependence on initial stellar mass and metallicity). For 
example, in the solar /super-solar metallicity range, beyond the Milky Way, only M31 has 
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been studied and just 48 WRs have been found in the whole galaxy f MasseylboOsI ). 

The Galactic center (GC) of the Milky Way (~8.0 kpc, iGhez et alj[2008[ ) provides an 
excellent laborator y for a systernatic s t udy o f EMSs with a apparent super-solar metallicity 
(~ 1.4 - 2x solar; llVIartins et allboOTl . |2008|) . Within the centr al 200 pc radius, the GC ac- 
counts for about 10% of t he ionizing photons (IFiger et al.ll2004j ) and ~25% of the entire WR 
population in the Galaxy (IMauerhan et al.ll2010d ). About half of the ioni zing radiation from 
the GC is attributable to the Arches, Quintuplet and Central clusters (ILang et al.l 120011). 
Thes e three massive star clusters have well-dete rmined ages (Arches: ~ 2 .5 Mvr. iFiger et al. 



2002 



Naiarro et al 



20041 : Martins et al.l I2008L Quintuplet: ^ 3.5 Mvr. iFiger et al. 



1999a: 



Liermann et al.ll2012l . and Central: ~ 6 Myr. lGenzel et al.ll2003l : iPaumard et al.ll2006[ ). Each 
has a total stellar mass of ~ 10^ Mq. They contain numerous identified EMSs, including 55 
WRs, at least two LBVs, and many OB supergiants. We can thus use these clusters to study 
the dependence of the evolution o f the EMSs on ini t ial ste llar mass and the impact that such 
stars have on their environment (ILang et al.l 119971 . 120011). In additiori , dozens of emission- 



line stars have been identified outside these clusters ( Cotera et al. 



Muno et al. 2006: Mikles et al. 2006: Mauerhan et al. 



2007, 



1999; 



20091 



i[ boiOah . 



Homeier et al. 2003 



While some of 

these stars were probably ejected from the clusters , others might form separately in small 
groups or even in isolation (IMauerhan et al.ll2010al ). We have recently carried out the first 
large-scale, high-resolution, narrow -band, Paschen-o; survey, using the NICMOS camera on 
the Hubble Space Telescope (HST) (IWang et al.ll2010[ ). In addition to the confirmation of the 
bulk of the known EMSs, the survey has led to the new detection of ~ 80 stars with signifi- 
cant excess emission in the 1.87 /xm band, which contains the HI Paschen-a and several Hell 
lines ( iDong et al.ll201ll ). These emission line stars should mostly be EMSs with strong stellar 
winds; MS stars and evolved low mass stars typically show absorption lines in their spec- 
tra. Indeed, follow-up spectroscopic observatio ns have already confirni ed 20 of the stars as 
EMSs such as WRs, LBVs and OB supergiants ( Mauerhan et al. 2010b ^); most of the other 
emission-line stars are fainter and are yet to be observed spectroscopically. By now, nearly 



100 W R. stars have been spectroscopically identified in the central 200 pc ( IMauerhan et al. 



2010d ). This s ample is dominated b y the late type WN/WC stars. In contrast, Westerlund 



1 (~ 4.5 Myr, ICrowther et al.l 120061 ) . another massive star formation region in the Galactic 
Disk, contains only 24 WR stars and half of the 16 WN stars are early-types (WN5-6). This 
difference is consistent with the higher metallicity of the GC than the Wd 1. Therefore, the 
detections and identifications of the massive stars in the GC now provide us with an excellent 
sample of EMSs or their candidates to study their individual and statistical properties and 
to infer relationships among the various types in a high metallicity environment. Further- 
more, EMSs are prominent signposts for recent massive star formation. From the study of 
the EMS populations in the GC, one can learn about the star formation mode, dynamics. 
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and history in an extreme Galactic nuclear environment, w hich is characterized by high gas 
temperature/density, strong magnetic field and tidal force (IMorris fc Serabyrull996l ). 



Here we report a multiwavelength analysis of EMSs in the GC. We describe the mul- 
tiwavelength data and the sample selection in § [21 including the new extinction correction 
and Paschen-a emission equivalent width (EW) calculation of the sample stars. In § [3l we 
present the near-IR, mid-IR and X-ray as well as emission line properties of the stars. In 
§ m we discuss the nature of the stars and their implications for the population, formation 
mode and history of massive stars. We summarize our results in § O Table 1 lists all the 
acronyms and abbreviations used in this paper. 



2. Sample Selection and Multiwavelength Data 



2.1. Sample Selection 

Table [2] lists 180 sources included in the present study; 152 of them are from a pri- 
mary list of the so-cal led Pascheii-a-em itting sources (PESs) uniformly identified in the 
HST/ NICMOS survey (IDong et al.ll201ll ). The identifications are based on the large 1.87 /im 



to 1.90 //m fiux ratios of these sources 



fpiSTN 
fpidON ' 



fpisTN and fpigoN are the intensities at 



1.87 and 1.90 /im, F187N, Paschen-a, on the Paschen-a line, F190N, on the adjacent con- 
tinuum). This ratio describes the Paschen-a excess for each individual star. For evolved, 
low-mass stars, which dominate the GC in a near-IR survey, this ratio is insensitive to 
the e xact stellar type (with < 1% variation) and should instead trace foreground extinc- 
tion (IDong et al.ll201ll ). The local ratios (f) towards different lines-of-sight are defined as 
the median average flux ratios of the closest 101 stars. All of the identifled PESs (ex- 
cept for some of the foreground stars, see section [3^ have 'signiflcant Paschen-a excess', 
Ng = ^-^ > 4.5, where atot includes both statistical and systematic photometric uncertain- 
ties of the ratio measurements. Ns is the signal-to-noise of the flux ratio of each source above 
the local ratio. Unfortunately, within the three clusters, this survey can be problematic in 
isolating individual stars, because of the limited spatial resolution of the NIC3 camera used 
for the survey. Its pixel size is 0.2", which undersamples the PSF (FWHM = 0.17" at 1.90 
/im). The large local background in the cluster cores further increases the photometric un- 
certainty (atot) and hence decreases the detection sensitivity, which may partly explain why 
more t han ten previously known WR st ars in the C e ntral cluster were not detected in our 
survey ( iDong et al .1120 111 ). Therefore, in iDong et al.l (120111 ). we also gave a secondary PES 
catalog, which includes the sources with 3.5<A^s<4.5. Compared to the primary catalog, this 
secondary catalog should include more spurious identiflcations. This catalog is only useful 
in the star clusters, in which we already know that there are many evolved massive stars. 
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since some of the massive stars would not be included i n the primary catalog because of the 
unusually high atot mentioned above (jPong et al.ll201ll ). 



Instead of using the secondary catalog to select the PESs in the three massive star 
clusters, we further refine the PES detection in the core regions of the Arches and Central 
clusters, utilizing archival HST/NIC2 F187N and F190N observations (Programs # 7250 
and 7364). By using the NIC2 data with its smaller pixel scale (0.074"), we can improve the 
photometric accuracy of the massive stars in these crowded regions. The basic parameters 
of these two sets of NIC2 observ ations are listed in Table El We reduce the data in the same 
way as was done in our survey ( iDong et al.ll201l[ ). In particular, we use the IDL program 
'Starfinder' f lDiolaiti et al.ll2000[ ) to detect the sources and to obtain their photometry with 
the P SF empirical l y extr acted from bright stars within the same observation. We use Eqn. 
(1) of iDong et al.l (1201 if ) to estimate the photometric uncertainty of the detected sources, 
including the Poisson fluctuations and the local background noise, as well as the uncertainty 
of the photometric calibration of individual filters. The undersampling uncertainty for NIC2 
data is neglected, because it is expected to be far less severe than for NIC3 data. With 
the improvements in the detection limit, in addition to updating the photom etries of the 



detec ted PESs (including the separation of one source , ID 122, in Table 3 of iDong et al. 



201ll into two separate stars, E48 and E51, as listed in lPaumard et al.ll2006l ). we find three 
and ten new PESs in th e Arch es and Central clusters, respectively (Table |2l), which are not 
reported by iDong et al.l ( 1201 ll . their Table 3). 



We also include 14 WC stars within our survey region that were identified by other 
spectroscopic observations, but are still not detected in our survey as significant PESs (MPl- 
MP14 in Table 121 Ng < 4.5; references for these stars are listed there). 



The completeness of our survey has been discussed in our previous paper (IDong et al. 



201ll ). Briefly, the typical 50% detection limit is shown to be about 17 mag, although it 
can be as high as 15.5 mag in the Central cluster region. With the NIC2 data used in the 
present work, the detection limit is further improved to 18.5 mag in the cluster regions. 
Because most of the 180 sources are significantly brighter than their local detection limits, 
we conclude that our EMS sample is nearly complete over the bulk of the survey field, except 
in regions with exceptionally high extinctions. 



2.2. Broad-band Near-Infrared Data 



Three broad-band near-IR data sets are used: 1) HST snapshot observations of the 
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Arches, Quintuplet and Central clusters (Programs 7250, 736 4, 7222 and 9457), 2) SIRIUS 



f lNagavama et all 120031 ). and 3) 2MASS fISkrutskie et all 120061 ) . The HST observations were 



targeted on the three clusters and their local background fields. Program 7364 utilized NIC2 
to map the Arches and Quintuplet clusters with three broad-band filters (FllOW, F160W 
and F205W): the covered fields are 1.4 arcmin^ (a total of 14 pointings) and 2.0 arcmin^ (20 
pointings), respectively. The Central cluster was mapped with one NIC2 (Program 7222) 
and four NIC3 pointings (Program7250) using F160W and F222M (covering a field of 2.8 
arcmin^). The medium filter (F222M) was used to avoid the saturation problem. There were 
also several HST/NIC2 observations with FllOW pointed toward a few diverse fields within 
the Central cluster from Programs 7222 and 9457. All of our PESs in the three clusters 
have counterparts in these broad-band observations. But they cover only about 1.5% of our 
HST/NICMOS Paschen-a GC survey field. 

We thus use the SIRIUS and 2MASS catalogues for those PESs not covered by the HST 
broad-band observations 0. SIRIUS covers our entire sample field in J (centered at 1.25 /^m), 
H (1.63 /im) and (2.14 /xm). For a few bright PESs, which are saturated in SIRIUS, we 
adopt the corresponding broad- band magnitudes (J: ce ntered at 1.24 ^m, H: 1.66 fim and K^: 



2.16 fj,m) provided by 2MASS (ISkrutskie et al.l 120061 ). We identify the closest broad-band 



source as the counterpart of each PES with a matching radius of 0.1". We also examine the 
individual matches by eye to make sure that the counterparts between the three catalogs are 
real matches. 

We convert the photometries of these catalogues into a common standard to facilitate a 
systematic analysis. Fig. [1] shows the transmission curves of the different systems in various 
filters. The filter response curves of the HST are much wider than those of the ground- 
based systems. We obtain approximate magnitude conversions using counterparts, or the 
closest sources in two different catalogs with offsets less than 0.2". For HST and SIRIUS, 
considering their much different filter systems, we first use a very crude match selection 
criterion: each source pair has a K-band magnitude difference less than 2 magnitudes. The 
H and filters of SIRIUS and 2MASS are very similar; we thus use a tighter criterion, 
requiring the magnitude difference to be less than 0.3. We then utilize least-squares fits 
to obtain the first-order color-dependent conversions to the SIRIUS magnitude system, and 
a 3 — cr clipping to remove the outliers. This 3 — a clipping effectively removes the SIRIUS 
sources with multi HST counterparts within the search radius, even in crowded regions such 



^Simultaneous 3-color InfraRed Imager for Unbiased Sur veys (SIRIUS) was tak en by the Infrared Survey 
Facility (IRSF) in South Africa, with a pixel scale of 0.45" (jNagavama et alJl2003l ). The survey includes the 
region |^| < 2 degree and |6| < 1 degree, with an angular resolution ^ 1.2" in the J band, better than that 
of 2MASS, -2". 
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as the Central cluster. Finally, we use least-squares fits again on the remaining source pairs 
to derive the following conversions: 



Jsu. = Jhst - 0.079(±0.006) (J^fST - Khst^w) 


-0.181(±0 


030) 


(1) 


Hs^r^ = HhST ' 0.011 (±0.015) (i///ST - KhST,F205w) 


- 0.228(±0 


024) 


(2) 


Ksiri = KhsT,F205W " 0.021 (±0.016) - KhsT,F205w) 


- 0.435(±0 


026) 


(3) 


Jsiri = Jhst - 0.159(±0.004)( J^,5t - Khst,F222m) 


- 0.247(±0 


031) 


(4) 


Hsiri = Hhst - 0.160(±0.004) (i:fH5T - KhST,F222m) 


- 0.246(±0 


Oil) 


(5) 


Ksiri = KhST,F222M ± 0.068(±0.003) (////^r — KhST,F222m) 


- 0.562(±0 


007) 


(6) 


Jsiri = J2mass ± 0.043(±0.001) ( J2mass ~ -^2mass) 


-0.021(±0 


003) 


(7) 


Hsiri = H2mass ~ 0.001 (±0 .001 ) (i^f2mass — K2mass) 


±0.075(±0 


002) 


(8) 


Ksiri = K2mass ± 0.070 (±0.001) (if2mass ~ -^2mass) 


- 0.075(±0 


002) 


(9) 



Fig. [2] compares the HST/ NICMOS, 2MASS and SIRIUS magnitudes before and after these 
conversions. The numbers of sources used to derive the equations above are also marked in 
each panel of Fig. [2l Table [2] includes the J, H and magnitudes of the PESs, all in the 
SIRIUS photometric system. 



2.3. Extinction Corrections 



Using the broad-band near-IR magnitude measurements, we improve our extinction 
correction for the PESs . In co ntrast to the use of the statistically constructed extinction map 
detailed in iDong et al.l ( 120111 ). here we directly estimate the extinction along the sightlines 
toward individual PESs, based on the J-H and/or H-K^ colors, if available. 

Table m lists the intrinsic colors, (J-H)o and (H-K3)o, adopted for the PESs. The large 
dispersion in the (H-Ks)o distribution of WR stars is mainly because their K^-band inten- 
sities can be substantially contaminated by the free-free emission from their strong stellar 
winds and/or the surrounding dust thermal emission. For those PESs with no available spec- 
troscopic identifications, we simply adopt the mean values obtained for O [(J-H)o =-0.11 and 
(H-K,)o =-0.1] and WR stars [(J-H)o=0.02±0.056 and (H-K,)o =0.19±0.182] (see Table. S]): 
(J-H)o =-0.045±0.121 and (H-K,)o = 0.045±0.327, the uncertainties of which include the 
color differences of O and WR stars. 



Following iNishiyama et al.l (|2006[ ). we adopt the relative extinction coefficients for the 
broad-band filters (J, H and K,) of SIRIUS toward the GC as Aj:Ah:Aks=1-0.573:0.331. 
Assuming a single power law for the extinction curve between H and K^, we then obtain the 
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slope, a=2.02, Afi9on=^-'^7^Aks, ^Fi87Af=l-306y4i^s and 

^Fi907v = 0.985(Aj - Ah) (10) 

Afi9on = 1-738{Ah-Aks) (11) 

^Fi907v = 90.56 xlog(^). (12) 

r 

The definition of f was given in § 12.11 

Fig. [3^ compares tlie Apig^j^ values obtained from Eqns. ITU] and [TT] for the PESs with 
available J, H and Kg measurements and known stellar types. The mean difference between 
^Fi90iv(J-H) and AFi90iv(H-K,) is 0.01±0.03 for OB supergiants, 0.02±0.04 for WN stars, 
and 0.16±0.09 for WC stars. For WC stars, AFi90Af(H-Ks) is overestimated because of the 
contamination by hot dust emission (§ 13. 4p in Kg. Therefore, we use only ApigQ]y{J-H.) 
to represent the foreground extinction of WC stars. For other PESs, we generally adopt 
the mean of ApigoNi'i-ii) and ApigoNiH-^s)- But for PESs with no J measurement, we 
adopt ApiQQ]\r{B.-Ks), if available; othe rwise, we just us e the extinction from Eq. [121 which 



represents the local average extinction (jPong et al.ll201l[ ). In Fig. Eb, we compare the Apig^j^ 



values derived from Eq. and Eq. [12] for the PESs with available J and H magnitudes, 
excluding the ten potential foreground stars identified in § 13.21 The similar extinctions 
derived from these two equations {ApiQQN{f)-ApiQQ]y{J-}i)=-0.1±0.07) suggest that most of 
these PESs are indeed within the GC. The extinctions so obtained are listed in Table [2 



2.4. Paschen-a Equivalent Widths 

We recalculate the Paschen-a equivalent widths (EWs) for those stars with improved 
magnitude and extinction measurements. The calculations follow the formula: 

??W fpiSJN ~ 1-015 * fpigoN x\ 

i.UiO * JpiQQN 

where fpis^j^f and fpigQj^ are the extinction-corrected fluxes of the PESs in the F187N and 
F190N bands, 6X ~0.191 A is the effective FWHM of the F187N fiher, and 1.015 represents a 
typical F187N to F190N continuum flux ratio, which is obtained for a KOIIl star and is insen- 



sitive to the exact stellar type, as long as there is no absorption lines in F187N (iDong et al. 



201l|). The EWs for the PESs are listed in Table [1 
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2.5. Mid-Infrared Data 

To further examine tlie dust emission from the PESs, particularly WC stars, we includ e 



the mid-IR measurements from the Spitzer/IRAC GALCEN survey (jStolovy et al.ll2006l ). 
We use a search radius of 1" to find the IRAC counterpart of each PES. Because of the 
limited spatial resolution of the IRAC camera, most of the EMSs within the cores of the 
Arches and Central clusters do not have mid-IR measurements. The mid-IR magnitudes of 
the PESs in the 3.6, 4.5, 5.4 and 8.0 /im bands are included in Table |2l The photometric 
uncertainties in the GALCEN source catalog are systematically underestimated by a factor 
of about 5 in the 3.6 and 4.5 jjm band^. Therefore, we correct for the uncertainties from 
the GALCEN catalog by this factor. 



2.6. Spectral Energy Distribution Analysis 



We conduct a simple spectral energy distribution (SED) fit to characterize the dust 
emission typically associated with WC stars. In our sample, 16 WC stars have available 
photometries in the six b ands: mj, uih, '>Pk„, " ^Fign/ y, [3.6] and [4.5]. We ad opt the ex- 
tinction laws provided by iNishiyama et al.l (120061 ) and llndebetouw et al.l ( l2005l ). to correct 



for the extinction in each of the six bands, scaled to ApigoN {Aj: Ah'- Aks'- ^3.6 tj.m'- 
2.38:1.36:0.79:0.43:0.34). Zero points of 835.6, 280.9, 179 .9 Jy are used to conver t 



A 



4.5 iMin 



the mag nitude into flux for F190N, [3.6] and [4.5] i As sugge sted bvlMauerhan et al.l (l2010al ) , 



because the photometry of SIRIUS and 2MASS are similar (INishiyama et al.ll2008l ). the zero 



points for the J, H and Kj, magnitudes of the 2MASS system are used for sources in the 
SIRIUS catalog (for which we find no zero point information). Five of these WC stars (P28, 
P34, P60, P94, PlOl) have relatively blue colors (K5-[3.6]<2) and do not show evidence for 
surrounding dust (those without 'd' in their type definition; Table [2]) . We first normalize the 
flux measurements of each of these WC stars to /j = 1 Jy and then median-average them to 
construct a "dust-free" SED template, which should include any potential free-free emission 
from the stellar winds as well as the emission emerging from the stellar atmospheres of these 
sources. The uncertainties of the median-average SED in the six bands are less than 20%. 
We fit the SED of the other 11 WC stars with this template, plus a diluted blackbody to 
characterize hot dust emission with an emissivity proportional to , /, 



"A4 exp(^)-l' 



m 



see 



http://www.astro.wisc.edu/sirtf/glin2_galcen_coinparison.pdf 



^littp://www.stsci.edu/list/nicmos/performance/pliotometry and Spitzcr/IRAC instrument handbook, 
Table 4.1 
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units of Jy, where a is the normahzation factor and h 
(see Fig. |4]). 



^ JCoudfrooii fc de JoneJEooi 



2.7. X-ray Measurements 



We utihze the measurements made by the deep Chandra ACIS survey (iMuno et al. 



20091 ) to search for X-ray counterparts of the PESs. The matching radius used is the smaller 
of 3(7^ (where is the X-ray source positional uncertainty) and 2" (which is used to reduce 
the source confusion). To estimate the probability for chance positional coincidences, we 
repeat the same matching procedure after each of eight systematic shifts of the X-ray source 
positions (i.e., the matches were actually done on a grid of 5RA C [-5", 0", 5"], 5Dec C [-5", 
0", 5"]). The mean number of chance positional coincidences in the eight shifts is found to 
be 7.5, which mostly occur in the core regions of the three clusters. Outside the clusters, 
the number of chance positional coincidences is only 1.8. 

Table [5] lis t s 35 P ESs with X-ray source counterparts. In addition to the source IDs listed 



m 



Muno et al.l (120091 ) . the table also includes X-ray source positional uncertainty (cr^), 0.5-8 



keV fluxes and hardness ratios (s ee the note to the tab le). P122a and P122b, members of the 
IRS 13 compact stellar complex (IMaillard et al.ll2004l ). are just 0.3" apart (compared to the 
Chandra on-axis FWHM ~0.5") and are listed to have the same X-ray source counterpart. 

Most of these PES/X-ray source matches should represent genuine counterparts. Table[5] 
includes the percent probability (P) for each PES/X-ray source match to be genuine. The 
probability was calculated by Mauerhan et al. (2009) using a match radius that is the 
quadrature sum of the positional un certainties of the Chandra and SIRIUS measurements 



Those matches studied previously (IMuno et al.l l2006l : iMikles et al.l l2006l : iMauerhan et al. 
2010al ) all have high probabilities (> 20 %), except for P147 (partly because of its large 
positional uncertainty: 0.9"). Because of saturation in the SIRIUS observation, the source 



P 38, which was id e ntified in our HST survey (IDong et al.l 1201 ll ). was not listed in Table 3 



of 



Mauerhan et al.l (120091 ) . and therefore does not have an existing 'P' value. iMauerhan et al. 



( 2OO9I ) also did not match any X-ray sources to the four PESs in the Central cluster (P123 



P125, P158 and P161 in Table [5]), because the SIRIUS catalog suf fers from confusion i n 
this crowded region. However, these sources are known WR stars (jPaumard et al.ll2006l ). 
Four other PES/X-ray sources (P105, P108, P133 and P139) have low 'P' (<15%) and have 
not yet been spectroscopically classified. Therefore, whether they are real matches remains 
uncertain (see also §|1]). 
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Results 



The locations of our 180 sample PESs are marked on the Paschen-a mosaic image of the 
GC in Fig. |5l A total of 30, 17 and 33 of the sources are located within 2rc of the Quintuplet, 
Arc hes and Cen t ral clus ters, respectively; the cluster radius, Tc, of each cluster is from Table 
5 of Figer et al. ( 1999a ) (Quintuplet: Ipc, Arches: 0.19 pc. Central: 0.23 pc) @. We refer to 
these sources as 'cluster' PESs and to the remaining 100 PESs as 'field' ones. 



3.1. Foreground Extinction Properties 



The extinction distributions of the PESs in different regions are depicted in Fig. |6l The 
means (standard deviations) of ApigoN are 3.32(0.56), 3.20(0.24), 3.71(0.80) and 2.98(0.89) 
for the Quintuplet, Arches, Central and 'field' PESs, respectively. The relatively large dis- 
persion for the Central cluster PESs is largely due to the exceptionally large extinctions 
of MP7 (E31, WC9) and MPIO (E58, WC5/6, with no available J-band measurement; Ta- 
ble [2]), most likely due to the extra extinction of their circumstellar dust. Excluding these 
two sources reduces the mean and standard deviation to 3.55 and 0.5. In general, the PESs 
within the individual clusters have very similar extinctions. The corresponding A^s (stan- 
dard deviation) for the Quintuplet, Arches and Central clusters are 2.61( 0.44), 2.52(0.19) 



and 2.79(0.39), which are consistent w ith existing me asurements: 3.1 (0.5 , iLiermann et al. 



20101), 2.13-4.14 (lEspinoza et al.ll2009[ ) and 2.54 (0.12. ISch5del et al.ll2010[ ). respectively 



Fig. [7] shows the J-H vs. H-K,; and vs. H-K^ diagrams. The color-color plot contains 
119 PESs with both J-H and H-K^ measurements, while the color- magnitude diagram is used 
for the 45 PESs with only H-K^ colors. Most of the PESs are located in the ranges of H- 
KsG [1, 3] and J-Hg [2, 4]. In Fig. [71 we also plot all of the stars in the SIRIUS catalog, which 
are mainly low mass giants, plus foreground MS stars. They should typically have similar 
intrinsic near-lR colors (on the Rayleigh- Jeans side of the SED). Therefore, the observed 
colors of such stars should trace the expected reddening path. In contrast, the PESs show 
colors that can be significantly redder than the path (e.g., up to 1-2 mag in H-K^), apparently 
a result of enhanced emission at long wavelengths, due to free-free and/or hot dust emission. 
The remaining 16 PESs do not have H or measurements. 



" ^The core radii of the Arches and Central clusters are consisten t with recent work 
pc lEspinoza et al.ll2009t) and Central: 0.22±0.04 p c JSchodel et al 



Arches cluster derived from dynamical simulations (jKim et al 



1999; 



Arches: 0.14±0.05 
20071 ). The tidal radiu s (^I pc) of the 



Portegies Zwart et al.l 120021 ) is larger. 



Using this latter radius, two additional PESs (P17 and P79) would be classified to be within the Arches 
cluster, which, however, would not qualitatively affect our results and conclusions. 



-VI- 



3.2. Identification of Foreground Stars 



The above extinction properties of the PESs can help us identify/confirm stars that are 
in the foreground of the GC. We find that nine stars have H -K.^ < 1 and J-H < 2 (if a vailable), 
including the two foreground OB supergiants suggested by lMauerhan et al.l (l2010d ). Table [6] 
includes the extinction-corrected F190N magnitudes and updated A^^ values (see § 12.11) of 
these foreground PESs. Four of them (P27, P102, P140 and P149) still have significant 
F187N excesses, with A^^ > 4.5. Although the A^^ values of the remaining five PESs are 
now < 4.5, they may still be foreground emission-line stars. For example, P38, with the 
second smallest revised A^^, s pectroscopically e xhibits several emission lines around 2.112- 
2.115 /im and is identified by iMauerhan et al.l (j2010d ) as an 04-61 star at a distance of 3.6 
kpc. 

The nine foreground PESs can be divided into two groups according to their extinction- 
corrected F190N magnitudes. The first consists of P38, P140 and P145, which have compa- 
rable foreground extinctions and intrinsic brightnesses. As both P38 and P140 are EMSs, we 
suggest that P145 is likely to be as well. The other six PESs form the second group. They are 
a factor of at least ~ 10 fainter and have comparable or smaller (but within a factor of ~ 3) 
extinctions, compared to those in the first class, indicating that they are typically clos er and 
substa ntial ly dimmer by rn o re tha n > 2.5 mag than P38, for example. According to iFiger 
(119951 ) and ICrowther et al.l (120061 ). even the dimmest WR (Mj^g ~-2.63) is no more than a 
factor of 10 fainter than the first group sources. Therefore, we suspect that the PESs in this 
second class are not EMSs (such as WR and OB supergiants) and are probably cata clysmic 
variables (CVs), which often show hydrogen emission lines ( iDhillon fc Marshlll995l ). CVs 
usually contain lo w-mass MS secondaries (< IMq) with M^^s in the 2-9 mag range (Fig. 9 of 
Hoard et al.ll2002l ). Therefore, with the extinction taken into account, only relatively nearby 
CVs (up to a distance of ~ 4 kpc) can be detected in our survey. 

When X-ray counterparts are available, their spectral properties can also be used to 
check the consistency of the identifications and even to provide clues to the nature of the 
sources. Table O summarizes the 35 matches between the PESs and known X-ray sources. 
Most of these matches have been examined with previous ground-based near-IR observations, 
and have been suggested to be massive coUiding-wind binaries within the GC (except for 
P38; iMauerhan et al.ll2010cl ). The X-ray har dness ratio HRO is ofte n used to distinguish 
foreground sources from those at the GC (e.g., Mauerhan et al. 2010a| j3). IMuuo et al. (2009) 
suggested that a source with HRO < -0.175 probably has a distance less than 4 kpc. There 
are four such EMSs, including P38, consistent with its foreground star identification. P145 
also has a soft spectrum, consistent with being a foreground massive star, according to its 
near-IR color above. The statistical uncertainties in the HRO values for P109 and P138 
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are too large to provide useful constraints. For example, the red near-IR color of P109 
{H — Kg = 1.6) indicates that it could be within the GC, instead of a foreground source. 

The nine stars with blue near-IR colors in Tabled plus one soft X-ray star (P138), are 
excluded from the discussion to follow, which focuses on the EMSs within the GC. 



3.3. Continuum and Paschen-a Emission Properties of GC EMSs 



We present the F190N magnitude distributions of our 170 GC EMSs in Fig. [HI The 
distribution becomes increasingly broad from the Arches (~2.5 Myr old). Quintuplet (~4 
Myr old), to Central (~6 Myr old) clusters. There are also systematic magnitude differ- 
ences among the distributions, which may be characterized by their mean F190N magnitude 
(standard deviation): ~-6.21(0.58), ~-6.93(2.09) and ~-5.26(1.33) for the Arches, Quintu- 
plet and Central clusters, respectively. The unusually large displacement of the Quintuplet 
cluster from the other distributions may be due to the presence of several very bright stars 
(AfpiQOAT < —8) due to substantial contamination by circumstellar dust emission: mainly one 
LBV and t he five Quintuplet proper members (QPM s, dusty WCs), which have unusually 
red colors ( iFiger et al.l Il999al : iLiermann et al.ll2009[ ). Excluding these six stars, the mean 
would decrease to ~-6.18(1.55). The standard deviation further decreases to 0.85, after the 
three dim stars with Mpi^Q^^ > —5 are removed. On the other hand, the stars in the Central 
cluster are fainter by ~ 1 -2 magnitudes. This offset is consisten t with the results derived 
from the Kg observations ( iFiger et al.ll2002l : iPaumard et al.l 120061 ) and is at least partly due 
to the larger age of the Central cluster. Therefore, the magnitude distribution shift in general 
seems to be linked to the age of a cluster. 

The magnitude distribution of 'field' sources spans a wide range and is possibly double- 
peaked. About 50% (46) of the 'field' sources have absolute 1.90 /im magnitudes similar to 
those in the three young clusters, while the other half are considerably fainter. Compared 
to those in the bright peak, the mean extinction of the PESs in the dimmer peak is smaller 
by ~ 0.4 magnitude, which cannot explain the large magnitude difference between these two 
peaks (> 2 mag). 

Fig. M presents an EW vs. fpigo^ plot for the PESs with identified stellar types in 
different regions. The stellar types are identified from previous ground-based spectroscopic 
observations (Table [2]). WNL stars in the Arches cluster (triangles) show an app arent linear 
relatio n between the EW and the stellar luminosity, as already demonstrated by iFiger et al. 
( 120021 ). They pointed out that this relation is expected from radiation-driven stellar winds. 
A similar linear relationship seems to hold for WNL stars in the Central cluster, excluding 
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the Ofpe/WN9 stars (filled squares). But these WNL stars are systematically dimmer by 
a factor of ~ 5 than those in the Arches cluster (see § 13.31 and Fig. |8]), implying that the 
EW is determined by not only the luminosity, but also the stellar age. The three WNL 
stars in the Quintuplet cluster are more similar to those in the Arches cluster, than to the 
Central cluster. In Fig. |9l most of 'field' WNL stars ('star' symbol) are located between the 
Arches and Central clusters. However, 'Ofpe/WN9' stars seem to be distributed differently 
and have substantially smaller EWs (< lOOA) than WNL stars of similar fpig^j^f- Both the 
Quintuplet and Central clusters are rich in WC stars, with a broad EW range. The three WC 
stars in the Quintuplet with EW > 250 A (P12, P13, P60, which are all WC8) have strong 
Hel (2.059 /xm) and Hell (2.19 /im) in their near-IR K^-band spectra flFiger et al.l Il999al : 



Liermann et al.ll2009[ ). But many WC stars have very weak line emission, which explains 



why 14 WC stars are not identified in our PES survey (see §|2]). Five of them are the QPMs, 
indicating that their stellar photospheres are completely blanketed by dust. Most of the OB 
supergiants also have very small EWs, indicative of their weaker stellar winds. 

In the three clusters, no EMS is detected to have fpig^j^f < 2 x 10^ jjJy and EW< lOOA 
at the same time. This is likely due to the relatively high flux detection limit in the cluster 
fields due to crowding. Fig. dU] shows the photometric uncertainty atot (see § 12. ip vs. /l-igoAr 
distribution of the detected sources within the three clusters, compared with the unclassified 
'field' EMSs, which have low fpigQ^ EWs (we use /^igoAr instead of fpidON to avoid the 
effects of the differential foreground extinction). Because of the high local stellar background, 
the photometric uncertainty in the Central cluster is much higher than those in the Arches 
and Quintuplet clusters. The photometric uncertainty for the Quintuplet cluster (which 
only has the low resolution of NICM0S/NIC3 F187N/F190N data) is also larger than that 
in the Arches cluster (for which NIC2 data are used; § 2.1). Therefore, the non-detection 
of the analogue of the field EMSs, which have relatively low brightness and EWs, in the 
Quintuplet and Central clusters could be due to their large atot- In the case of the young 
Arches cluster, however, massive stars are still mostly on the MS, except for a few very 
massive ones (> 60Mq) that may have just turned into bright WNL stars or OB supergiants. 



3.4. Dust Emission in the Mid-IR 

We present a color-color diagram for the 36 EMSs with available H, and [8.0] mag- 
nitude measurements in Fig. [TTJ Limited by the angular resolution of Spitzer/IRAC (~2"), 
these PESs are all in the 'field', although a few of them are located at the peripheries of 
the Quintuplet and Arches clusters. Most of the WN stars and OB supergiants are only 
slightly redder in Ks-[8.0] than average field stars. The three OB supergiants (P35, PlOO, 
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P112) that have unusually red color {Ks — [8.0] > 5) are all associated with local, strong, ex- 
tended Paschen-a features, so their [8.0] flux density is presumably dominated by Polycyclic 
Aromatic Hydrocarbon (PAH) emission. While all WN stars are constrained to [2,4], a sig- 
niflcant number of WC stars have very red K5-[8.0] color. Six of the nine WC stars with red 



K,^-[8] ( >3) in Fig. [TT] have been suggested to be associated with hot dust (IMauerhan et al. 



2010al ]d: iLiermann et al.ll2009l : iFiger et al.lll999al ). 



Fig. [H shows an EW vs. K,-[3.6] plot for 68 EMSs. We can see that WNL stars/OB 
supergiants and WC stars seem to follow different trends. For WNL stars, the EW increases 
with Ks-[3.6]. But, the trend is opposite for WC stars, at least for those with low EWs. 
The EWs of the OB supergiants show a very weak, possibly positive dependence on K^- 
[3.6]. These different trends probably indicate different relative contributions from various 
sources of emission. The free- free emission from the strong winds of th e WNL stars/OB 
supergiants should dominate in the mid-infrared (jWright fc Barlowlll975l ). resulting in the 
positive EW vs. Ks-[3.6] correlation. In contrast, the strong dust emission from WC stars 
may dilute the Paschen-a line emission and dominate the mid-infrared. For example, QPMs 
have very red infrared spectra and appear nearly featureless. Indeed, dusty spirals have 
been resolved spati ally around two of these WC stars based on high-resolution adaptive 
optics observations (iTuthill et al.ll2006l ). The enhancement of the hot dust emission at the 
longer wavelengths suppresses the line contrasts. As a result, a signiflcant fraction of WC 
stars may thus not be detected as PESs. 

The results from the SED flts to the 11 WC stars (§ 12. 6p with enhanced mid-IR emission 
are presented in Table [7] and Fig. |4l The combination of the "dust-free" template and 
the diluted dust blackbody flts the bulk of the WC stars well. For MPl, we notice that 
the Spitzer/IRAC source catalog includes an additional source that is only 1" away and is 
detected only at 4.5 /im, but not at 3.6 /xm. This source may thus represent a spurious 
detection and may have caused an underestimate of the intensity of MPl at 4.5 /im in 
the Spitzer/IRAC catalog. Therefore, we have excluded the 4.5 /im data point from the 
SED fltting of MPl in Table [7] and Fig. HI Both the table and flgure show that the dust 
contribution at 1.90 //m varies signiflcantly from one source to another and reaches more 
than 30% for the QPMs. After removing the dust contribution from the F187N and F190N 
flux densities, we recalculate the EWs, the results of which are also listed in Table [71 For 
the WC stars not detected in our survey (MP*), their EWs are increased, though still very 
small (< 30A). Since dus ty WC st ars are presurnably formed within colliding- wind binary 
systems (ICrowtherl 120071 . see also iTuthill et al.l l2006l ). the presence of massive (OB MS) 
companions with flat spectra and even Paschen-a absorption lines could also account for 
the low EWs. The fltted dust temperatures are in the ran ge 550-1000 K. We der i ve th e 
dust mass in these WC stars (see Table. [7]) using Eqn. 2 of ICoudfrooij fc de Jongl (119951 ): 
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Md=5.1 X 10-^i5^D2A^(ei-^^^i°*/^''^'' - 1) Mq, where (in units of mJy) and D (Mpc) are 
the specific blackbody flux at frequency u and the distance to the source, respectively. In 
order to exclude the contamination from the stellar atmosphere at short wavelengths, we 
use Si, at 4.5 /im, except for MPl (at 3.6 //m). The dust mass in each of the five QPMs 
(MPl, MP2, MPS, MPS and MP6) is usually large, 10"^ - 10"^ Mq, while the other six WC 
stars have dust masses of order 10~^ — 10~^ Mp^, similar t o the WC stars studied in the solar 
neighbourhood (jVeen et al.lll998l : iMarchenko et al.lll999l ). 



4. Discussion 

In this paper, we construct a catalog of 180 EMSs in the Galactic center region, in- 
cluding many identified WN, WC, LBV and OB supergiants. Based o n the brightness and 



EW distributions of such stars, as discussed in lMauerhan et all (j2010d ). this catalog should 
include nearly all of the WN stars in our survey area, while a significant fraction of faint WC 
stars may still be missing. Even in the core regions of the Arches and Central clusters, while 
the Paschen-a survey itself may be quite incomplete because of the crowding and nebulosity, 
the dedicated NIC2 observations of higher spatial resolution and ground-based spectroscopic 
observations should have identified the bulk (if not all) of the WN stars. 

Based on the above results, we can gain some insight into the nature of those still 
unclassified PESs, especially those outside the clusters. Fig. [13] compares the PESs with and 
without spectroscopic classifications in the EW vs. f°QQ plot. For ease of the discussion, 
we empirically divide the plot into three parts, each of which is populated mainly by one 
single stellar type: the WNL-dominated region contains 67% of the WNL stars, the OB 
supergiant-dominated region contains 96% of the OB supergiants, and the WC-dominated 
region contains 42% of the WC stars. 53% of the WC stars with very small EWs fall into the 
'OB supergiant' region. In Fig. [121 we find that the WC stars with red color (larger Ks-[3.6]) 
tend to have smaller EWs. Therefore, we propose that the SEDs of WC stars falling into 
the 'OB supergiant' region could be substantially contaminated by the thermal emission of 
hot circumstellar dust, which may affect the flux even at 1.90 fim. 

The unclassified PESs have either low EWs and/or low intrinsic brightnesses. 54 out of 
67 such PESs are within the 'OB supergiant' region (Fig. [T3|) . One in the Central cluster 
with EW=343A falls within a tight cluster of WNE stars in Fig. [121 and is detected in th e 
HST/NIC2 observations, but not in the HST/NIC3 Paschen-a Survey JPong et al.lboilh . 
because of the confusion limit in the very crowded region. The other 12 unclassified PESs 
are all in the 'WC region. 
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The unclassified PESs are scattered around in the plots of Fig. [TT] and Fig. [T2j Their 
MIR properties suggest that they include sources from various subgroups, such as dusty WC 
stars and OB supergiants. For example, P44, P115 and P141 are very red [Kg — [8.0] > 4 
and Kg — [3.6] > 3) and are unrelated to any local diffuse emission in the Paschen-a map. 
Therefore, they are most hkely dusty WC stars. Several other unclassified PESs fall into 
the 'OB supergiant' region. Spectroscopic identification is needed to determine their exact 
stellar types. 

The six previously unclassified PES/X-ray source matches (Table [S]) are all located 
outside the three clusters. They have 0.5-8 keV fluxes of ~ 4x 10~^ photons cm~^ s~^, dimmer 
than most of the other X-ray sources in Table [51 Because of the low signal-to-noise ratios of 
the sources, we only use their X-ray colors for spectral characterization. In § 13. 2^ HRO values 
are used to check the consistency of classifying P138 and P145 as foreground stars. Fig. [T^ 
further presents the relationship between the HR2 and the photon fluxes of these X-ray 
counterparts. Unlike HRO, HR2 is a measure of the spectral charac teristics at relatively hig h 
energies, and is less affected by foreground extinction. As stated in lMauerhan et al.l (j2010al ). 
the X-ray counterparts of the EMSs are mostly softer and brighter than the majority of the 
X-ray point sources toward the GC, which are typically CVs and have relatively hard spectra, 
with thermal temperature > 2 keV. Among the six new detections, P105 and P133 have HR2 
(< 0.4) similar to the other known EMSs with X-ray counterparts, suggesting that they may 
have a similar origin. On the other hand, P108 and P139, with relatively hard HR2 (> 0.5), 
fall into the CV region in Fig. [HJ However, as we i nentioned in ^ 13. 2 [ after correcting for 
the distance modulus and the extinction (^^'=2.22, iDong et al.l 1201 if ), the CVs in the GC 
should be dimmer than 18.7 magnitude, which is well below the source detection limit of our 
survey (see Fig. 14. of iDong et al.ll201ll ). Therefore, P108 and P139 are probably unrelated 
to the X-ray point sources, consistent with their poor matching probability (Table |5]). 



The number of 'fie ld' EMSs is comparable to that within the three clusters (IMauerhan 



20081 : iDong et al.l 1201 ll ). Fig. |8] shows that the absolute magnitude distribution of these 



'field' EMSs is much wider than those of individual clusters and shows two peaks, perhaps 
suggesting two distinct episodes of star formation. The absolute magnitudes of more than 
half of the EMSs fall between the peaks of the Arches and Central clusters. Some of these 
EMSs are WNLs and are mostly located between the loci defined by WNLs within the Arches 
and Central clusters in the EW vs. F190N plot of Fig. |9l These facts suggest that these 
EMSs have ages between those of the Arches and Central clusters. The remaining fainter 
unclassified 'field' PESs are presumably older. Some of them have magnitudes comparable 
to those PESs in the Central cluster that are WC stars {M^s G [—3, —5]). Therefore, some 
of these unclassified 'field' PESs could also be WCs and trace the stellar population with 
ages comparable to that of the Central cluster. 
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Although some fraction of the field EMSs may represent the EMSs that have been ejected 
from the three clusters via dynamical processes, at least some, especially those young ones, 
likely formed in small groups or even in isolation. Probably the most convincing evidence 
for this sce nario is the appar ent associations of a few of EMSs with local HII regions, such 
as HI, H2 ( IZhao et all Il993l ). and several others in the Galactic West of Sgr A* (Fig. |5l 
Fig. [15]). In addition, the unusual red colors (K5-[8.0]) of the three OB supergiants also 
indicate their associations with nearby molecular clouds (§ 13. 4p . Such associations should be 
mostly real; the possibility for chance coincidence of the EMSs with dens e clouds is smal l, 
because the clouds have a small volume filling factor (~ 1%) in the GC flOka et al.ll2005l ). 
The relatively small size (<25", i.e. <1 pc, see Fig. [T5|) and faintness of the associated HII 
regions further indicate that they are ionized by small groups of massive stars or probably 
by the EMSs alone. The EMSs in HI, H2 and H5 are spectroscopically identified as OB 
supergiants. The total extinction-corrected Paschen-a luminosity of the HI and H2 regions 
(see Fig. [T5b) is 6.1x10^^ ergs s ~^, corresponding to a total ionizing photon rate of ~ 
4.1 X 10^^ (jScoville et al.ll200ll ). This rate can be accounted for by the two enclosed OB 
supergiants (P35 and P114; Qiyc ~ 10^^'^ each), consistent with the formation scenario 
that some of the field PESs formed in small groups or in isolation. 



5. Summary 

In this paper, we have investigated the properties of 180 PESs by using multi- wavelength 
observations from near-IR to X-ray. These sources have been selected using their significant 
intensity enhancements in the F187N narrow band, based on our HST/ NICMOS survey of 
the GC and previous HST snapshot observations of the Arches and Central clusters. The 
multi-wavelength data analysis has enabled us to probe the overall spectral characteristics 
of individual sources, including a direct estimation (hence correction) of line-of-sight extinc- 
tions/absorptions. Our main results follows: 

• 10 sources are identified to be foreground stars. While a few of them are apparently 
EMSs, the others with intrinsically low luminosities could be CVs. 

• For WNL stars, the Paschen-a line EW is well correlated with the intrinsic magnitude 
and color Ks-[3.6], which traces the free-free emission, suggesting that the winds of the 
WNL are radiation-driven. However, the EW dependence on the magnitude appears 
to be sensitive to the stellar age. 



• The above correlation seems to hold for WC stars in which the Paschen-a line EW is 
greater than about 100 A. But those WC stars having lower EWs show anti-correlations 
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with the magnitude and the Ks-[3.6] color, which can be explained by an increased 
continuum level due to the strong hot dust emission from the stellar winds. 

The presence of dusty stellar winds is clearly shown in the IR SEDs of a subset of the 
WC sample, showing intensity upturns or peaks in the mid-IR. From the modeling of 
the SEDs, we estimate that the dust has a characteristic temperature of ~ 750 K and 
that each star has a dust mass in the range from 10~^ Mq to a few x 10~^ Mq (see § l3.4p . 
The upper end of this range is much higher t han had previously been known for WC 



stars ( IVeen et alll998l : lMarchenko et alll999l ). which may have important implications 



for understanding the production of dust, e specially in very young starburst galaxies 



at high redshifts ( iDwek fc Cherchneflal201ll ). 



We show that PESs may be roughly typed according to their positions in the EW vs. 
magnitude diagram. This typing works 'best' for WN stars. Both EW and magni- 
tude ranges of WC stars are broad. Some WC stars have very low EWs, apparently 
due to enhanced continuum emission from hot circumstellar dust. They cannot be 
distinguished from OB supergiants based only on the EWs and magnitudes. We have 
tentatively assigned those unclassified PESs as dusty WC/OB supergiants. Our detec- 
tion of WR stars with relatively strong Paa line emission should be quite complete. But 
it is possible that many Paa-faint WRs, especially WCs associated with substantial 
amounts of dust, remain to be detected. 

Nearly half of the PESs are located outside the three clusters. The magnitude distribu- 
tion of these 'field' sources shows two distinct peaks. One of the two peaks, containing 
roughly half of the sources, has a mean magnitude similar to those of the EMSs in the 
three clusters, while the other is substantially dimmer. The unclassified PESs in this 
latter peak have F190N magnitudes similar to or fainter than those of the WCs in the 
Central cluster. This suggests that the dimmer peak traces a stellar population with 
an average age similar to, or older than, that of the Central cluster. 

Considering that the volume filling factor of dense clouds in the GC is very small, we 
suggest that a few of the field PESs are physically associated with compact HII regions 
(with sizes < 25" or 1 pc) and formed in situ, indicating the operation of a mode of 
massive star formation in small groups or even in isolation. 
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Fig. 1. — Filter transmission curves of various instruments. Notice that the filters of 
HST/NIC3 are much wider than those of SIRIUS and 2MASS, and that the J filter of 
SIRIUS is narrower than that of 2MASS. 



Table 1:: Acronym/ Abbreviation List 



CV 
EMS 



Cataclysmic variable 
Evolved massive star 
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Table 1:: Acronym/ Abbreviation List 



EW 


Equivalent width 


FWHM 

_L V V J- J-i-V J. 


l^'nll wiHtVi at brilf rnayimmn 


GC 


Galactic center 


HST 


Hubble SDace TelescoDe 


IRAC 


Infrared array camera 


ISM 


Interstellar medium 


LBV 


Luminous blue variable 


MS 


Main sequence 


NICMOS 


Near infrared camera and multi-object spectrometer 


PAH 


Polycyclic aromatic hydrocarbon 


PES 


Paschen-CK emitting source 


QPM 


Quintuplet-proper member 


SED 


Spectral energy distribution 


SN 


Supernova 


WR 


Wolf-Rayet 


WN 


Nitrogen sequence Wolf-Rayet 


WNE 


Early type Nitrogen sequence Wolf-Rayet 


WNL 


Late type Nitrogen sequence Wolf-Rayet 


wc 


Carbon sequence Wolf-Rayet 


WCE 


Early type Carbon sequence Wolf-Rayet 


WCL 


Late type Carbon sequence Wolf-Rayet 



Table 2. Catalog of Paschen-a emission-line stars 



Name R.A. 


Dec. 








Magnitude 
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EW 


Ns 


Counterpart 


Type 


Location 


(J2000) 


(J2000) 
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H 


K 


[1.87] [1.9] 


[3.6] 


[4.5] 


[5.4] 


[8.0] 















PI 

P2 
P3 
P4 
P5 
P6 
P7 
P8 
P9 
PIO 
Pll 
P12 
P13 
P14 
P15 
P16 
P17 
P18 
P19 
P20 
P21 
P22 
P23 
P24 
P25 
P26 
P27 
P28 
P29 
P30 
P31 
P32 
P33 
P34 
P35 
P36 
P37 
P38 
P39 
P40 
P41 
P42 
P43 
P44 
P45 
P46 
P47 
P48 
P49 



266.62478 
266.59921 
266.51501 
266.55606 
266.55408 
266.56301 
266.56643 
266.56293 
266.56312 
266.56304 
266.56896 
266.56323 
266.56316 
266.58167 
266.51338 
266.47853 
266.46028 
266.45712 
266.45255 
266.45865 
266.45895 
266.47249 
266.54169 
266.50251 
266.49295 
266.49541 
266.48141 
266.49067 
266.53998 
266.53377 
266.52612 
266.37803 
266.33953 
266.46061 
266.36926 
266.38126 
266.44762 
266.40831 
266.34458 
266.35029 
266.40788 
266.38541 
266.30814 
266.25555 
266.25250 
266.28706 
266.31271 
266.34449 
266.34120 



-28.78001 
-28.80299 
-28.78606 
-28.81641 
-28.82014 
-28.82693 
-28.82714 
-28.82480 
-28.82568 
-28.82631 
-28.82547 
-28.82821 
-28.82761 
-28.83606 
-28.81622 
-28.78699 
-28.82543 
-28.82372 
-28.82840 
-28.82393 
-28.82411 
-28.82693 
-28.92566 
-28.90761 
-28.87223 
-28.89392 
-28.90196 
-28.91267 
-28.95375 
-28.97338 
-28.98747 
-28.87674 
-28.86082 
-28.95726 
-28.93473 
-28.95466 
-29.04884 
-29.02624 
-28.97893 
-29.01633 
-29.10817 
-29.08277 
-29.07728 
-29.04208 
-29.10650 
-29.11563 
-29.15230 
-29.18235 
-29.19988 



13.4(0.02) 
16.3(0.05) 
17.5(0.11) 
14.6(0.02) 
17.5(0.07)'' 
13.5(0.06)'' 
14.2(0.06)'' 
14.1(0.06)'' 
14.3(0.06)'' 
14.5(0.06)'' 
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15.0(0.02) 



16.3(0.03) 
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14.5(0.02) 
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Fig. 2.— The difference between the SIRIUS and HST/NICMOS or 2MASS magnitudes as 
a function of the near-IR colors before (blue) and after (red) the magnitude transformation 
for the HST/NICMOS magnitudes and 2MASS. The values in each lower right corners are 
the number of sources used to determine the magnitude transformations. 
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Fig. 3. — (a) Comparison between the 1.90 /im extinctions of the PESs having available J, 
H and measurements, derived separately using J-H and H-K^ colors. The black solid line 
represents Ai;'i90Ar(J-H)=A^i90Ar(H-Ks). (b) Comparison betweeh the A^igoA? values of the 
PESs with available J and H magnitudes, derived separately using J-H and f (see § 12.11 The 
ten potential foreground PESs identified in § 13.21 which have blue color: J — H < 2, are 
excluded). The black solid line represents Airi90Ar(J-H)=Ai?i90Ar(f). 
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Fig. 4.— The SED fitting result for 11 WC stars in six bands: J, H, F190N, K^, 3.6 fim and 
4.5 fim. The dotted fines represent tfie "dust-free" template of WC stars, wfiile tfie dasfied 
lines characterize the diluted dust blackbody contributions. The solid lines are tfie sum of 
tfie two components. QPMs fiave been flagged as sucfi. All of tfiem fiave strong mid-IR 
emission. 
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Fig. 5. — 180 PESs overlaid on the mosaic image of Paschen-O! intensity ('diamond': WNL 
stars, 'square': WC stars, 'triangle': WNE stars, 'circle': OB supergiants, 'star symbols': 
LBV stars). The 'pentagon' symbols are the stars without available spectroscopic identifi- 
cations. The 'white' symbols indicate the sources with X-ray counterparts (see § 12. 7p . 
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Fig. 6. — The extinction distributions of the PESs in the Quintuplet (Yellow), Arches (Blue), 
Center (Red) and 'field' (Green). PESs with ApiQQ^ < 1.7, i.e., H-Ks<l (see Eqn. 11) are 
considered to be foreground stars. 
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Fig. 7. — Left panel: Infrared color-color diagram of 119 PESs with available J, H and 
magnitude compared to the field stars from the SIRIUS survey (small black dots). Right 
panel: Infrared magnitude-color diagram of 45 PESs without J detections, compared to the 
field stars from the SIRIUS survey (small black dots). Thick black lines delineate the region 
occupied by foreground stars (H-K <1 and J-H < 2) 
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Fig. 8. — F190 magnitude distributions before (left) and after (right) the extinction and 
distance corrections (see § I3.3p . The histograms for Arches, Central and 'field' are shifted 
upward by 10, 20 and 30 counts, respectively. 
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Fig. 9. — Paschen-a equivalent width (EW) vs. extinction-corrected F190N flux densities, 
fpidoN^ plots of stars of type WNL (a), WC (b), WNE (c) and OB supergiant (d). In (a), 
the fllled symbols are the Ofpe/WN9 stars. The solid lines in (b)-(d) represent the positive 
correlation between the EW and fpigoN, seen in (a) for the WNL stars in the Arches cluster. 
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Fig. 10. — Comparison of the uncertainty of the flux ratio f^H™ )^ ^toti (see § l2.ip between 
the normal sources within the three clusters and the unclassified PESs, which have low EW 

and /l-igoAT- 
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Fig. 11.— E-Ks vs. K,-[8.0] plot for 36 PESs and field stars from the Spitzer/IRAC GAL- 
CEN survey (black dots). 
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Fig. 12. — EW vs. Ks-[3.6] plot for various types of PESs. 
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Fig. 13. — The combined Paschen-a EW vs. /^iqotv P^^t of stars of various types. Lines are 
drawn to roughly show the regions dominated by WNL, WC and OB supergiants. The stars 
with larger symbols are the ones with X-ray counterparts. 
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Fig. 14. — HR2 vs. Photon flux plot f or the 34 X-ray cou nterparts of the PESs. The blue dots 
are from the Chandra X-ray Catalog (IMuno et al.ll2009l ). Most of them have high HR2 (>0) 
and low photon flux (< 10~^ cm~^ s~^) and should be CVs. The 'star' symbol represents 
the counterparts that have been reported in the literature, while the 'circle' symbols are the 
new X-ray counterparts. We also label the positions for different types of spectra , either 
power-law with index P or thermal plasma at temperature T, as in lMauerhan et al.l ( l2010a ) 
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Fig. 15. — HST and Spitzer images of four HII regions in F190N (left panel), Paschen-o; 
(middle) and 8.0 ^m (right). From top to bottom: (a) H5, (b) HI and H2, (c) (l,b) = (- 
0.06,0.02), (d) (l,b) = (-0.13,0.0) 
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3.0 


112.2(3.6) 


18.9 


FA_9,Maul0a_12 


WN8-9h 


A 


P87 


266 


45954 


-28 


82136 


15.7(0.07)'^ 


12.5(0.06)'' 


10.6(0.06)'' 


10.9(0.01) 


11 


6(0 


01) 




















3.1 


185.2(4.5) 


26.1 


FA_12,Blu01.25 


WN7-8h 


A 


P88 


266 


46121 


-28 


82284 


15.5(0.06)'^ 


12.5(0.05)'' 


10.7(0.06)'' 


11.1(0.01) 


11 


7(0 


01) 




















2.9 


155.1(4.1) 


16.7 


FA_14,Blu01.12 


WN8-9I1 


A 


P89 


266 


46151 


-28 


82118 


15.3(0.06)'^ 


12.3(0.05)'^ 


10.7(0.06)'' 


11.4(0.01) 


11 


6(0 


01) 




















3.0 


47.1(2.8) 


8.1 


FA.15,Blu01.8 


04-6If 


A 


P90 


266 


46057 


-28 


82231 


15.7(0.06)'^ 


12.6(0.05)'^ 


10.8(0.06)'' 


11.3(0.01) 


11 


8(0 


01) 




















3.0 


116.8(3.7) 


13.8 


FA_16,Blu01.19 


WN8-9h 


A 


P91 


266 


48072 


-28 


85726 


15.4(0.19) 


12.5(0.04) 


11.0(0.07) 


10.5(0.03) 


11 


5(0 


02) 


9.4(0 


05) 


8 


7(0 


05) 


8 


5(0 


03) 




2.5 


319.9(17.1) 


18.8 


M07.1,Maul0a_16 


WN5-6b 


F 


P92 


266 


52345 


-28 


85881 


12.7(0.04)* 


9.3(0.03)* 


7.5(0.03)* 


7.1(0.02) 


7.4(0.02) 






5 


8(0 


03) 










3.3 


79.4(6.4) 


11.9 


MaulOb 


LBV 


F 


P93 


266 


51671 


-28 


89814 






15.1(0.08)* 


15.3(0.03) 


15 


6(0.04) 




















2.8'' 


61.2(11.1) 


5.5 






F 


P94 


266 


51080 


-28 


90388 


16.4(0.04) 


13.5(0.02) 


11.6(0.02) 


11.5(0.02) 


12.5(0.02) 


9.9(0 


07) 


9 


3(0 


07) 


9 


1(0 


03) 




2.9 


287.1(13.2) 


21.8 


Maul0c_14 


WC9 


F 


P95 


266 


49765 


-28 


88075 


14.9(0.02) 


12.1(0.01) 


10.5(0.02) 


11.1(0.02) 


11.2(0.02) 


9.6(0 


06) 


9 


2(0 


07) 








9.8(0.22) 


2.7 


24.8(5.1) 


4.8 






F 


P96 


266 


45238 


-28 


83491 


17.6(0.10) 


13.3(0.02) 


11.0(0.02) 


11.0(0.02) 


11 


8(0.02) 


9.1(0 


04) 


8 


5(0 


04) 


8 


1(0 


02) 


7.7(0.06) 


4.0 


229.2(10.0) 


22.8 


Cot.l 


Ofpe/WN9 


F 


P97 


266 


44891 


-28 


84692 


15.3(0.19)* 


12.3(0.09)* 


10.7(0.06)* 


11.3(0.02) 


11.5(0.02) 


9.5(0 


06) 


9 


1(0 


06) 


8 


9(0 


03) 




2.9 


47.8(6.6) 


7.6 






F 


P98 


266 


42197 


-28 


86325 


14.8(0.02) 


11.6(0.01) 


9.9(0.01) 


10.4(0.03) 


10.6(0.02) 


8.6(0 


04) 


8 


0(0 


04) 


7 


8(0 


02) 


7.6(0.05) 


3.2 


52.6(8.2) 


6.6 


Maul0c_9 


04-6If-l- 


F 
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Table 2 — Continued 



Name R.A. 


Dec. 








Magnitude 










A 


EW 


Ns 


Counterpart 


Type 


Location 


(J2000) 


(J2000) 


J 


H 


K 


[1.87] [1.9] 


[3.6] 


[4.5] 


[5.4] 


[8.0] 















P99 


266 


41118 


-28 


860^8 






1 3 Ql'n 041 


1 4 '^^n 091 


1 1 (n 091 


PlOO 


266 


42634 


-28 


87976 


14.7(0.02) 


11.7(0.01) 


10.1(0.03) 


10.5(0.03) 


10.9(0.02) 


PlOl 


266 


42704 


-28 


88140 


16.3(0.04) 


13.2(0.01) 


11.2(0.02) 


11.3(0.02) 


12.1(0.02) 


P102 


266 


43349 


-28 


88800 


15.0(0.02) 


13.3(0.03) 


12.4(0.04) 


12.5(0.02) 


12.8(0.02) 


P103 


266 


50689 


-28 


92091 


13.4(0.01) 


10.7(0.01) 


9.1(0.05) 


9.4(0.02) 


9.7(0.02) 


P104 


266 


53042 


-28 


95485 








15.9(0.05) 


16.3(0.07) 


P105 


266 


47070 


-28 


92692 




14.9(0.04) 


13.4(0.10) 


13.8(0.02) 


14.0(0.02) 


P106 


266 


46776 


-28 


94615 




16.1(0.08) 


14.5(0.13) 


14.4(0.02) 


14.7(0.02) 


P107 


266 


41391 


-28 


88919 


14.7(0.02) 


11.8(0.01) 


10.2(0.01) 


10.5(0.02) 


10.9(0.02) 


P108 


266 


44602 


-28 


94612 




15.5(0.07) 


14.0(0.09) 


14.4(0.03) 


14.7(0.03) 


P109 


266 


46088 


-28 


98877 


15.4(0.02) 


12.5(0.01) 


10.9(0.01) 


10.8(0.02) 


11.6(0.02) 


PllO 


266 


45393 


-28 


98382 




15.7(0.04) 


14.3(0.03) 


14.8(0.03) 


15.0(0.03) 


Pill 


266 


40061 


-28 


94403 


15.6(0.03) 


12.3(0.01) 


10.4(0.02) 


10.4(0.02) 


11.1(0.02) 


P112 


266 


40753 


-28 


95450 




i O. i(U.Uz ) 


iU. ( {yj.yjZ ) 


ll.b(U.U^) 


TO c\( {\ 
LZ.\J\\}.\)Z ) 






jyuoo 








±41. Z(^U. Uo ) 








P114 


266 


'^86^8 


-28 


93794 


1 4 7^0 09"^ 






11.2(0.02) 


11 iTi noi 


P115 


266 


32590 


-28 


89079 


16.4(0.06) 


13.8(0.04) 


11.9(0.02) 


12.7(0.02) 


12.9(0.02) 


P116 


266 


29639 


-28 


91955 


13.7(0.02) 


13.0(0.01) 


12.8(0.01) 


12.7(0.02) 


12.8(0.02) 


P117 


266 


41477 


-29 


00973 




12.8(0.03)'' 


10.6(0.02)'' 


10.8(0.01) 


12.1(0.00) 


P118 


266 


41443 


-29 


00881 




1 fifTi riO\^ 


11 ^f'nnil'^ 

J-J-.0(U.UJ- ) 


19 i^n nil 


19 sfn nil 


P119 


266 


41377 


-29 


00853 




14.4(0.03 ) 


Iz.z(O.Oz) 


1^.D(U.U1 ) 


lo.o(O.Ul ) 


P120 


266 


41410 


-29 


UU9^9 




14.8(0.04 ) 


lz.e(O.Oo ) 


lo.o(U.Ul ) 


lo.9(U.Ul ) 


P121 


266 


41777 


-29 


00756 


14.0(0.06) 


10.9(0.02)'' 


8.9(0.01)'' 


10.0(0.01) 


10.1(0.00) 


P122a 


266 


41584 


-29 


00829 


16.0(0.07)'* 


12.8(0.03)'' 


10.5(0.02)'' 


11.3(0.02) 


12.0(0.01) 


P122b 


266 


41583 


-29 


00821 




13.8(0.04)'' 


11.1(0.02)'' 


12.1(0.04) 


12.9(0.02) 


P123 


266 


41728 


-29 


00458 


16.7(0.06)'* 


13.6(0.03)'' 


11.4(0.01)'' 


11.6(0.02) 


12.4(0.03) 


P124 


266 


41635 


-29 


00504 


16.8(0.06)'' 


13.8(0.02)'' 


11.8(0.01)'' 


12.3(0.01) 


12.9(0.01) 


P125 


266 


41555 


-29 


00739 


19.6(0.15)'' 


14.7(0.04)'' 


11.5(0.02)'' 


13.7(0.02) 


13.9(0.01) 


P126 


266 


41560 


-29 


00646 


17.5(0.06)'' 


14.2(0.02)'' 


12.2(0.01)'' 


12.8(0.01) 


13.3(0.01) 


P127 


266 


41714 


-29 


00621 


17.6(0.08)'' 


14.2(0.05)'' 


11.4(0.02)'' 


12.7(0.04) 


13.0(0.04) 


P128 


266 


41608 


-29 


00616 


18.9(0.09)'' 


15.5(0.04)'' 


13.2(0.02)'' 


13.9(0.02) 


14.5(0.02) 


P129 


266 


41717 


-29 


00767 


15.1(0.06)'' 


11.7(0.02)'' 


9.7(0.01)'' 


10.8(0.01) 


11.1(0.00) 


P130 


266 


41777 


-29 


00814 


18.1(0.11)'' 


14.2(0.04)'' 


11.7(0.02)'' 


12.3(0.01) 


13.3(0.02) 


P131 


266 


41984 


-29 


00772 


17.9(0.11)'' 


14.7(0.03)'' 


12.8(0.01)'' 


13.2(0.01) 


13.8(0.01) 


P132 


266 


41774 


-29 


00936 


17.5(0.07)'' 


14.3(0.03)'' 


12.0(0.02)'' 


12.5(0.02) 


13.2(0.01) 


P133 


266 


40359 


-29 


02154 






13.9(0.08)* 


15.1(0.05) 


16.0(0.10) 


P134 


266 


31975 


-28 


97363 


16.7(0.06) 


13.1(0.03) 


11.1(0.02) 


11.4(0.03) 


12.1(0.02) 


P135 


266 


32402 


-28 


97220 




16.1(0.10) 


14.1(0.06) 


14.8(0.03) 


15.1(0.03) 


P136 


266 


36759 


-29 


05754 




14.8(0.03) 


12.7(0.02) 


13.2(0.02) 


13.3(0.02) 


P137 


266 


31749 


-29 


05437 


11.6(0.03)* 


9.2(0.03)* 


7.9(0.03)* 


7.7(0.04) 


8.3(0.03) 


P138 


266 


33870 


-29 


13177 








14.3(0.02) 


14.5(0.02) 


P139 


266 


31898 


-29 


09621 




15.4(0.05) 


13.5(0.03) 


14.0(0.02) 


14.2(0.02) 


P140 


266 


24782 


-29 


09057 


8.8(0.03)* 


7.7(0.04)* 


7.0(0.03)* 


6.9(0.02) 


7.2(0.02) 


P141 


266 


28870 


-29 


13783 




15.3(0.05) 


11.7(0.02) 


13.0(0.02) 


13.2(0.02) 


P142 


266 


33105 


-29 


17609 




15.9(0.06) 


14.8(0.13) 


14.8(0.03) 


15.0(0.03) 


P143 


266 


34964 


-29 


17353 








15.8(0.04) 


16.2(0.04) 


P144 


266 


30924 


-29 


19493 








15.3(0.03) 


15.6(0.03) 


P145 


266 


32051 


-29 


20505 


10.0(0.04)* 


8.7(0.05)' 


7.9(0.03)* 


8.1(0.02) 


8.2(0.02) 


P146 


266 


22471 


-29 


09574 








13.8(0.03) 


14.6(0.04) 



8.5(0.04) 
9.4(0.05) 



7.7(0.07) 
8.7(0.05) 



6.5(0.01) 
7.6(0.03) 7.1(0.04) 6.7(0.01) 



8.9(0.05) 

9.6(0.05) 

8.8(0.04) 
7.8(0.04) 



5.6(0.04) 



8.5(0.07) 

9.0(0.05) 

8.3(0.04) 
6.9(0.03) 



8.1(0.03) 

8.6(0.02) 

8.0(0.02) 
6.1(0.01) 



4.6(0.01) 
5.6(0.01) 

7.9(0.04) 

7.8(0.06) 
5.1(0.01) 



i. 9(0. 04) 8.1(0.04) 7.1(0.01) 6.2(0.01) 



9.4(0.05) 8.8(0.05) 
6.3(0.03) 



11.0(0.11) 
6.8(0.03) 



6.5(0.02) 
7.6(0.04) 



6.1(0.02) 
6.4(0.03) 



7.3(0.03) 
8.7(0.03) 

11.0(0.14) 
6.1(0.01) 



5.9(0.01) 
5.5(0.01) 



5.8(0.03) 
7.4(0.03) 

9.3(0.15) 
5.8(0.02) 



6.0(0.01) 
5.0(0.01) 



7.4(0.03) 7.3(0.03) 7.3(0.01) 7.1(0.02) 



4.2 
3.0 
3.0 
1.6 
2.9 

3.0'' 
2.7 
2.7 
2.9 
2.5 
2.8 
2.3 
3.1 
4.3 
3.4 
2.6 
2.9 
0.5 
3.5 
3.4 
3.7 
2.7 
3.1 
3.4 
4.2 
3.3 
3.1 
5.0 
3.2 
3.3 
3.4 
3.3 
3.9 
3.1 
3.4 

3.7'' 
3.4 
3.5 
3.7 
2.2 

2.8'' 
3.2 
1.3 
6.3 
1.8 

3.3'' 

3.0'' 
1.3 

3.0'' 



200.7(11.4) 
77.5(9.0) 
201.4(9.3) 
40.9(5.6) 
56.2(5.8) 

104.7(23.2) 
37.6(6.1) 
62.9(7.9) 
83.0(7.5) 
49.0(8.3) 

195.3(10.6) 
34.9(8.3) 
171.9(8.6) 
110.0(7.4) 
37.5(5.5) 
40.6(5.5) 
24.9(6.0) 
13.4(4.9) 

415.6(3.8) 
172.5(2.8) 
268.4(4.1) 
158.4(4.1) 
39.4(1.5) 
201.0(6.2) 

241.7(15.3) 

218.1(13.2) 
152.5(2.7) 
67.5(5.2) 
114.4(2.3) 

63.3(13.3) 
144.1(8.1) 
64.0(1.8) 
306.2(9.3) 
162.2(4.7) 
209.5(6.5) 

238.6(43.6) 

174.0(12.3) 
70.0(9.3) 
33.0(5.5) 

135.3(15.4) 
41.4(6.5) 
43.0(6.2) 
58.3(5.9) 
57.7(6.0) 
37.0(7.8) 
74.1(14.8) 
77.1(10.9) 
16.3(4.9) 

224.5(20.0) 



17.4 
8.2 
21.0 
8.8 
9.7 
4.5 
6.2 
8.1 
11.4 
6.4 
19.2 
4.7 
20.0 
13.8 
6.9 
7.1 
5.1 
6.0 
58.5 
23.9 
39.7 
22.9 
7.6 
25.2 
15.1 
16.5 
16.8 
9.6 
13.6 
4.7 
12.5 
13.6 
28.1 
23.9 
14.3 
5.5 
13.8 
7.3 
5.8 
8.7 
6.4 
6.6 
11.3 
6.7 
5.6 
5.0 
7.1 
4.7 
11.2 



MaulOc.lO 
Hoi 



Cot_5 
Cot_2,MaulOa_15 



Mik06.01,Maul0a_9 
Cot_6 



E79 
E74 
E81 
E82 
E39 
E51 
E48 
E88 
E83 
E56 
E66 
E71 
E68 
E20 
E40 
E78 
E65 



Muno06.01,Maul0a_3 







F 


04-6If+ 




p 


WC8-9 




F 






F 


OI 




F 






F 






F 






F 


B[e] 




F 






F 


WN7 




F 






F 


WN8-9h 




F 


B[e] 




F 








04-6I 










p 






p 


Wipe/ w In y 






WN8 




C 


WN7 




c 


WC8/9 




c 


Ofpe/WN9 




c 


WN8 




c 


WC9 




c 


WN8/9 




c 


WN8/WC9 




c 


Ofpe/WN9 




c 


WN8 




c 


WC8/9 ? 




c 


WC9 




c 


Ofpe/WN9 




c 


WN5/6 




c 


WC9 




c 


WN8 




c 






F 


WN7-8h 




F 






F 






F 


Ofpe/WN9 




F 






F 






F 


B0I-B2I 




F 






F 






P 






F 






F 






F 






F 
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Table 2 — Continued 



R.A. Dec. Magnitude A EW Ns Counterpart Type 

(J2000) (J2000) J H K [1.87] [1.9] [3.6] [4.5] [5.4] [8.0] 







-29 




lo. l(^U.Uo ) 


11 IZI (\ A1\ 
i.Z.K3\\J.\Jd. ) 


il.l(^U.U4 ) 


10.8(^0.0^ } 


1 1 |5 fr\ AO \ 
ll.Ol^U.Uz ) 








2^ 


T3n Q (■ 1 A 1\ 

zoii.oyLV.z ) 




iviauiua_z 


WN7 




P148 


zoo .OViiZV 








±0 . yv.i-V ) 


1/1 '7( c\ ^Q^ 


10 . Z\V .Vo) 


±0.0 (^U.Uo ) 


















P 


P149 


266.31115 


-29 


25314 


Id 1 (C\ C\^') 
L^. L y \j .\J 1. } 






12.4(0.02) 










1.1 




15.0 






F 


P150 


266.31271 


-29 


24339 




14.5(0.02) 


12.7(0.03) 


13.1(0.02) 


13.7(0.02) 


11.4(0.12) 


10.8(0.11) 


10.0(0.05) 


3.1 


143.7(8.5) 


16.5 






F 


P151 


266.29083 


-29 


23696 


15.6(0.03) 


13.0(0.01) 


11.1(0.02) 


11.6(0.03) 


12.1(0.02) 


9.4(0.05) 


8.7(0.05) 


8.4(0.03) 


8.0(0.04) 2.5 


105.6(10.0) 


10.6 


Maul0c_4 


WC9?d 


F 


P152 


266.24459 


-29 


25164 








14.9(0.04) 


15.1(0.03) 








2.8*' 


57.9(12.3) 


4.7 






F 


P153 


266.45859 


-28 


82313 


16.7(0.07)'' 


13.2(0.06)'' 


11.3(0.06)'' 


12.1(0.01) 


12.2(0.01) 








3.4 


25.6(2.6) 


4.6 


FA_13,Blu01.31 




A 


P154 


266.46092 


-28 


82199 


16.3(0.07)'' 


13.2(0.05)'' 


11.5(0.06)'' 


12.2(0.01) 


12.5(0.01) 








3.0 


48.1(3.1) 


5.7 


FA.27,Blu01.16 




A 


P155 


266.45930 


-28 


82120 


16.6(0.06)'' 


13.4(0.05)'' 


11.6(0.06)'' 


12.5(0.01) 


12.6(0.01) 








3.3 


28.3(2.8) 


5.2 


FA.22,Blu01.27 


04-6I 


A 


P156 


266.41716 


-29 


00808 


14.9(0.06)'' 


12.0(0.02)'' 


10.0(0.01)'' 


10.9(0.01) 


11.3(0.01) 








3.1 


77.1(2.2) 


16.0 


E23 


Ofpe/WN9 


C 


P157 


266.41684 


-29 


00747 


15.0(0.06)'' 


12.2(0.02)'' 


10.1(0.01)'' 


11.1(0.01) 


11.2(0.00) 








3.1 


35.5(1.5) 


10.5 


E19 


Ofpe/WN9 


C 


P158 


266.41704 


-29 


00868 


15.8(0.07)'' 


12.1(0.02)'' 


10.0(0.01)'' 


11.1(0.01) 


11.3(0.00) 








3.5 


40.0(1.8) 


6.6 


E41 


Ofpe/WN9 


C 


P159 


266.41545 


-29 


00827 


18.0(0.12)'' 


14.4(0.03)'' 


12.2(0.02)'' 


13.4(0.04) 


13.9(0.02) 








3.6 


113.3(11.9) 


9.1 


E60 


WN7? 


C 


P160 


266.41757 


-29 


00789 




14.7(0.04)'' 


12.7(0.02)'' 


13.6(0.03) 


13.9(0.03) 








3.6 


82.9(11.0) 


7.1 


E27 


O9.5-B0.5 I 


c 


P161 


266.41823 


-29 


00643 


18.8(0.09)'' 


15.0(0.03)'' 


12.7(0.02)'' 


13.4(0.03) 


14.0(0.01) 








3.8 


154.8(10.4) 


12.9 


E70 


Ofpe/WN9 


C 


P162 


266.41644 


-29 


00608 








13.3(0.04) 


14.3(0.05) 








4.6*' 


342.7(28.4) 


11.5 






C 


P163 


266.41820 


-29 


01005 




15.3(0.05)'' 


13.1(0.07)'' 


13.7(0.04) 


14.2(0.02) 








3.4 


121.2(11.2) 


6.9 


E76 


WC9 


C 


P164 


266.41566 


-29 


00701 


19.1(0.10)'' 


15.5(0.04)'' 


13.1(0.02)'' 


14.2(0.02) 


14.5(0.02) 








3.8 


71.3(6.6) 


9.3 


E61 


WN7 


C 


P165 


266.41777 


-29 


00685 


18.6(0.10)'' 


15.5(0.05)'' 


13.1(0.02)'' 


14.1(0.03) 


14.6(0.03) 








3.1 


114.2(11.0) 


9.5 


E59 


WC9 4^ 


C 


MPl 


266.56614 


-28 


82928 


15.2(0.10)' 


10.5(0.05)' 


7.4(0.04)* 


8.7(0.02) 


8.6(0.02) 








4.6 


2.2(4.6) 


-1.5 


Lie.l9,FQ_211,Maul 


wc8/9d -i-oe^ 


Q 


MP2 


266.56131 


-28 


82795 


13.6(0.08)'' 


9.4(0.07)'' 


6.4(0.08)'' 


7.8(0.02) 


7.7(0.02) 








4.2 


-3.0(4.4) 


-2.2 


Lie.42,FQ_231,Maul 


WC9d + OB| 


Q 


MPS 


266.55892 


-28 


82683 


15.7(0.08)'' 


11.1(0.08)'' 


7.9(0.08)'' 


9.5(0.02) 


9.4(0.02) 








4.5 


-1.7(4.5) 


-2.3 


Lio_75,FQ.243 


WC9?d 


Q 


MP4 


266.56418 


-28 


82624 


15.9(0.09)'' 


12.1(0.07)'' 


9.5(0.07)'' 


10.5(0.02) 


10.5(0.02) 








3.8 


12.7(4.8) 


1.8 


Lio_79,FQ.250 


WC9d 


Q 


MP5 


266.56167 


-28 


82616 


14.9(0.08)'' 


10.8(0.07)'' 


7.8(0.08)'' 


9.3(0.02) 


9.2(0.02) 








4.1 


0.5(4.5) 


-1.5 


Lie_84,FQ.251 


WC9d 


Q 


MP6 


266.55970 


-28 


82537 


17.9(0.10)'' 


13.2(0.08)'' 


9.4(0.09)'' 


10.8(0.02) 


10.7(0.02) 








4.7 


-1.1(4.5) 


-2.1 


Lie.l02,FQ_258 


WC9?d 


Q 


MPT 


266.41633 


-29 


00742 


19.3(0.24)'' 


13.2(0.03)'' 


9.9(0.02)'' 


12.2(0.01) 


12.1(0.01) 








6.0 


7.0(1.7) 


0.4 


E31 


WC9 


c 


MP8 


266.41742 


-29 


00814 


17.7(0.13)'' 


13.6(0.03)'' 


10.8(0.02)'' 


12.5(0.02) 


12.5(0.01) 








4.0 


11.7(3.3) 


2.6 


E32 


WC8/9 


c 


MP9 


266.41654 


-29 


00724 


19.0(0.15)'' 


14.6(0.04)'' 


11.7(0.02)'' 


13.5(0.02) 


13.5(0.01) 








4.3 


16.7(3.8) 


3.4 


E35 


WC8/9 


C 


MPIO 


266.41609 


-29 


00676 




15.3(0.05)'' 


11.0(0.02)'' 


13.7(0.01) 


13.6(0.01) 








6.4 


-3.1(2.2) 


-1.2 


E58 


WC5/6 


C 


MPll 


266.41898 


-29 


00796 


18.9(0.11)'' 


15.7(0.06)'' 


13.6(0.02)'' 


14.9(0.14) 


14.7(0.04) 








3.2 


-25.8(21.6) 


-1.0 


E72 


WC9? 


c 


MP12 


266.41862 


-29 


01009 




13.8(0.04)'' 


11.6(0.03)'' 


12.4(0.09) 


12.7(0.03) 








3.5 


88.7(25.5) 


3.5 


E80 


WC9 


c 


MP13 


266.32978 


-29 


05614 




13.7(0.02) 


10.4(0.01) 


10.7(0.02) 


10.8(0.02) 








5.3 


28.6(5.2) 


2.9 


Maul0a_5 


WC9d 


F 


MP 14 


266.57427 


-28 


83530 


14.2(0.06)* 


10.4(0.06)* 


7.9(0.04)* 


9.2(0.02) 


9.2(0.02) 








3.8 


11.2(4.8) 


1.3 


Maul0a_23 


WC9d 


F 



Note. — R.A. and Dccl. arc in degrees. J, H and K band magnitudes are mainly from the SIRIUS catalog jNishiyama et al.|f2b06l ). except for those with the superscript 'h' and 't', which are from 
the HS T/NIC2 snap s hot o bservations and 2MASS catalog (§ |2.2t . The extinction A values at l.QO^tm are from the near-IR color s, except for those with the ^superscript 'b', which are from the ratio 
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1 iDong et all j201ll1 (§ [2731. EW 

I 20091 ) f Lie.*). iPaumard et al.l l2006|)fE*') [bT 

I 2007l HMau07-*). iMauerhan et al.l l20l'0al 'i f MaulO; 



in units of A. The 
um et al.l <200l|') f BluOl.*) 
[ Mauerh an 



pcct roscopically identifie d count erpart names and typ es are from ^^iger e^^Ll il999al ) (F Q,*) , [F iger et 
aO Jl99gl 'i(' CQt_*1.lHomeier et al.l 120031' ) f Ho*) . iMuno et al.l <200a 'l f Muno06-*) ■ iMiklcs et al.i l| 



et al, 



Cotera et ; 



iLiermann 


et 


al. 


iMauerhan 


et 


al. 



■ <2010bl' )f Maul0b_ .. ^ „ ^ , ^ 

reference, except for P6 and P7, which are from iFige^^^^LLT^^QQElT r instead of iLiermann e^^i!Lr ' l2009l 7. The 'Location' column divides the sources into four groups: those inside the three clusters {'Q' - 
Quintuplet, 'A'- Arches, 'C - Center) and outside ('F'; §[^. The MPl-14 are not selected out by our method (N^ < 4.5, see ^ |2. H . Instead, they are WC stars identified by previous ground-based spectroscopic 
observations. We include these WC stars into our catalog to complete the sample of the evolved massive stars in the GC 
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Table 3. Previous HST NIC2 observations 



Program ID 


Observation Data 


Filter 


exposure 


7250 


Aug 31st, 1999 


F187N 


1215 sec 


7250 


Aug 31st, 1999 


F190N 


1215 sec 


7364 


Sep 14th, 1999 


F187N 


1024 sec 


7364 


Sep 14th, 1999 


F190N 


1024 sec 



Table 4. The intrinsic near-IR colors of massive stars 



Typo 


Description 


(J - H)o 


(J - H)q Error 


(H - K,)o 


{H ~ Ks)a Error 


WN5-6 weak 


WN5, WN6 


0.03 


0.012 


0.16 


0.012 


WN4-7 strong 


WN4b, WN5b, WN6b 


0.10 


0.006 


0.27 


0.006 


WN7-9 weak 


WN7(h), WN8(h), WN9(h) 


0.02 


0.052 


0.11 


0.053 


WNlO-11 


Ofpo/WNQ, WNIO, WNll 


0.06 


0.049 


0.13 


0.11 


WC5-7 


WC5, WC6, WC7 


0.04 


0.029 


0.58 


0.012 


WC8 


WC8 


0.05 





0.38 





WC9 


WC9 


-0.03 





0.26 





LBV 


LBV 


-0.07 





-0.05 





OB supcrgiant 


O If+, B[cl 


-0.11 





-0.1 






Note. — The colors of the WNs and WCs arc from ICrowther et aLl 1200^). According; to ^ iger et al.| 
j 19981 ). we assume the colors of LBV arc similar to 09.5I <Koo^inee3|l983l ) . iMartins Sz Plej i2006l ) provides 
the colors of OB supergiant 
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Table 5. X-ray counterparts of the Paschen-a emitting sources 



X-ray ID 



(arcscc) X-ray flux 



P (%) 



P9 


Lie_96,Maul0a_19 


CXOGC174615. 


1- 


■284932 





.51 


1 


,4c- 


-06 


PIS 


FA^.MaulOa.ll 


CXOGC174549. 


,7- 


■284925 





.41 


1 


,6e- 


-06 


P35 


Cot_4,Maul0a_7 


CXOGC174528. 


,6- 


■285605 





.48 


9, 


,0e- 


-07 


P36 


Maul0c_6 


CXOGC174531. 


,4- 


■285716 





.50 


2 


,4e- 


-07 


P38 


Maul0c_7 


CXOGC174537. 


,9- 


■290134 


0, 


.33 


1 


,4c- 


-06 


P39 


Maul0a_6 


CXOGC174522. 


,6- 


■285844 


0, 


.63 


2 


,3c- 


-07 


P75 


M07.2,MaulOa_22 


CXOGC174617. 


,0- 


■285131 


0, 


.56 


2 


,4c- 


-06 


P80 


FA_6,Maul0a_13 


CXOGC174550. 


,4- 


■284922 





.33 


4, 


,0c- 


-05 


P83 


FA_7,MaulOa_14 


CXOGC174550. 


,4- 


■ 284919 


0, 


.33 


2 


,3c- 


-05 


P86 


FA_9,Maul0a.l2 


CXOGC174550. 


2- 


■ 284911 


0, 


.33 


3, 


,0c- 


-05 


P91 


M07.1,Maul0a_16 


CXOGC174555. 


,3- 


■285126 





.49 


1 


,7c- 


-06 


P105 




CXOGC174552. 


9- 


■285537 





.59 


4, 


.7e- 


-07 


P108 




CXOGC174546, 


,9- 


■285645 





.56 


3, 


.6c- 


-07 


P109 


Cot.2,Maul0a_15 


CXOGC174550. 


,6- 


■285919 


0, 


.43 


2 


Ic- 


-07 


Pill 


Mik06.01,Maul0a.9 


CXOGC174536. 


1- 


■285638 





.37 


2 


,0c- 


-05 


P114 


Maul0a_8 


CXOGC174532 


7- 


■285617 





.50 


7, 


,2c- 


-07 


P122a 


E51 


CXOGC174539, 


7- 


■290029 





.30 


1 


,1c- 


-05 


P122b 


E48 


CXOGC174539, 


7- 


■290029 





.30 


1 


,1c- 


-05 


P123 


E88 


CXOGC174540, 


1- 


■290016 





.40 


2 


.4c- 


-07 


P126 


E56 


CXOGC174539. 


,7- 


■290026 


0, 


.30 


9, 


,8c- 


-06 


P131 


E7S 


CXOGC174540. 


,7- 


■290027 


0, 


.35 


4, 


Ic- 


-07 


P133 




CXOGC174536. 


,8- 


■290117 





.46 


1 


,6c- 


-07 


P134 


MaulOa_4 


CXOGC174516, 


7- 


■285824 





.65 


3, 


,6c- 


-07 


P137 


Muno06_01,Maul0a_3 


CXOGC174516, 


1- 


■290315 





.44 


3, 


,4e- 


-06 


P138 




CXOGC174521. 


,3- 


■290754 





.75 


7, 


,1c- 


-07 


P139 




CXOGC174S16, 


,6- 


■290545 


3 


.75 


4, 


.9e- 


-07 


P145 




CXOGC174516, 


,9- 


■291217 


1 


.32 


1 


.2c- 


-06 


P147 


Maul0a_2 


CXOGC174508, 


,9- 


■291218 





.87 


3, 


,4c- 


-06 


P1S8 


E41 


CXOGC174540, 


,0- 


■290030 





.30 


1 


,4c- 


-05 


P161 


E70 


CXOGC174540. 


,4- 


■290024 





.38 


2 


,1c- 


-06 


MPl 


Lie_19,FQ.211,Maul 


CXOGC174615. 


,8- 


■284945 





.49 


1 


,8c- 


-06 


MP2 


Lio_42,FQ.231,Maul 


CXOGC174614. 


,6- 


■284940 





.62 


8, 


,9c- 


-07 


MP12 


E80 


CXOGC174540, 


4- 


■290036 





.35 


9, 


.le- 


-07 


MP13 


Maul0a_5 


CXOGC174519. 


1- 


■290321 


0, 


.51 


1 


,0c- 


-06 


MP 14 


Maul0a_23 


CXOGC174617. 


7- 


■285007 


0, 


.50 


1 


,6c- 


-06 
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1.00 
1.00 
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1.00 
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-10.32 
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-10.37 
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9.00 
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-0.163 
0.561 
-0.123 
-O.53I 
-O.57I 
O.31I 
-0.053 
O.O2I 
-0.043 

o.oel 
-0.00^ 



9.00 
+0.57 

-0.48 

+0.95 
9.00 

+0.05 
0.07 



9.00 
+0.65 
0.55 
+ 0.94 
9.00 
+ 0.02 
0.02 
+ 0.19 
0.18 
+ 0.03 
0.03 
+0.03 
0.03 
-9.00 
-0.91 
0.04 
0-04 
-9-00 
-9^00 
+ 1^07 
9^00 
+ 0^36 
0.44 
+ 0.07 
0.07 
+0.76 
0.93 
13+0.75 
"■^■^-1.05 
+0.37 
0.81 
+0.39 
0.42 
+0.04 
0.06 



-0.21 
1.00 
0.07 

-0.03 
0.19 
0.19 

-1.00 
0.461 

-9.00 
1.00 
0.05 

-0.29 

-0.25 



-0.63 
0.54 
0.40 



„ (,,+0.16 



0.57 
-0.36 
0.21 
0.34 
0.06 



+ 0.36 
0.35 
+ 0.32 
0.31 
+ 0.18 
0.18 

+ 0.26 
0.26 
+ 0.44 
0.34 
+ 0.26 
0.25 
+ 0.07 
0.07 
+ 0.09 
0.09 
+ 0.08 
0.08 
+ 0.28 
0.25 
14+0-31 



-0.51 
0.11 
-0.32 



-0.073 
-O.33I 
-O.27I 
-O.lll 
-O.27I 
-O.O2I 
-0.64+ 



0.28 
+ 0.30 
0.34 
+ 0.58 
0.46 

20+°-°^ 
"■^"-0.02 
+ 0.19 
0.21 
„ 2_+0.03 
"■^'-0.03 
„ 2-+0.03 



0.59 
-0.14 



-0.13'^ 



-0.03 
-9.00 
-9-00 
5+0.04 
-0.04 
1+0.81 
-9.00 
- + 0.60 
' -0.67 
, + 0.44 
-0.40 
, + 0.10 
-0-10 
■9.00 
-1^20 
0.53+°-'>° 



-9.003 
-0.061 

i.ool 

O.25I 
-0.263 
-O.33I 
-1.00" 



-1.00 
0.37 
0.38 



0-42 
-9^00 
-1.34 
+ 0.24 
0.24 
+ 0.03 
0.03 



+0.26 
0.25 
+0.33 
0.35 
+0.29 
0.32 
+0.24 
0.24 
+0.25 
0.26 



-0 04+°-15 
"-"^-0.16 



-0.22 
-1.00 
0.09 
0.16 
-0.48 



+0.23 
0.23 
9.00 
0.82 
+ 0.20 
0.25 
+ 0.16 
0.17 
+ 0.27 
0.28 



43.5 
31.8 
41.4 
44.8 
N/A 
21.4 
31.9 
84.0 
73.5 
70.9 
41.0 
13.0 
N/A 
-1 
70.6 
39.8 
54.1 
54.1 
N/A 
N/A 
11.4 

9.2 
31.4 
79.8 
57.6 

0.4 
-1 
10.3 
N/A 
N/A 
25.9 
26.8 
27.7 
28.1 
48.1 



Note. — The X-ray ID, flux (photons s"^ cm^^ , 0.5-8 keV) and X-ray colors (HR) arc from iMuno et al.l i2009l ). HR=-j|^. 
For HRO, HRl and HR2, h (s) represents the X-ray photons in the 2.0-3.3 keV (0.5-2.0 keV), 3.3-4.7 keV (2.0-3.3 keV), 4.7-8.0 keV 
(3.3-4.7 keV). '1' ('-1') means that the source is not been detected in the soft (high) band. '-9' means that the source is not detected 
in both h and s bands, cr^ is the X-ray source positional uncertainty. The 'P (%)' column lists the percent probability for the 
infrared object as real counterpart of the X-ray corresponding source as listed in Table 3 of lMauerhan et al ] j2009t) . 'N/A' means 
that the PES is outside their adopted match radii of the X-ray sources. means that the PES has a red near-IR color (H-Ks>0.9), 

while the X-ray source is soft (HRO<-0. 175) . We fi nd that P145 actually has a blue near-IR color {H-Ks^O.8), based on 2MASS 
measurements, although a red color was adopted in iMauerhan et al ] l2009tl . using the SIRIUS catalog, in which the J and H band 
fluxes are saturated. 
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Table 6. Modified A^^ for tlie foreground stars 



Name Counterpart Extinction-free "^iJ^iooiV ^F190JV Revised Ns 



First Class 


P38 


Maul0c_7 


7.6 


1.5 


2.9 


P140 


MaulOc.l 


5.9 


1.3 


9.8 


P145 




7.0 


1.3 


3.3 






Second Class 






PI 




12.2 


0.5 


3.4 


P27 




13.8 


1.5 


7.0 


P58 




10.2 


1.6 


4.1 


P102 




11.2 


1.6 


7.2 


P116 




12.4 


0.5 


2.7 


P149 




11.8 


1.1 


13.6 



Note. — Ng is the parameter which we use to identify the PESs in ^ 12. 1| 
In iDong et al ] i2011l ). we only considered the sources with Ng > 4.5 as 
PESs. We divide the nine sources with blue near-IR colors into two class 
according to their extinction-corrected free itt-fwon (sec ^ |3.2| . The first 
class should be foreground massive stars, while the second class could be 
foreground CVs instead. 



Table 7. Dust properties in the WC stars 



Name Tdust (K) F190N dust percentage Unc EW Corr EW Dust Mass (Mq) 



P63 


834±53 


15.5± 4.7 


99.5±6.9 


120. 5± 11.4 


3.5 


xlO-^ 


-8 


P66 


574±91 


2.2± 2.2 


206.4±13.4 


211. 9± 14.8 


1.2 


xlO" 


-7 


P73 


727±44 


8.1± 2.9 


49.7±8.1 


56.1± 9.5 


6.6 


xlO" 


-8 


P151 


944±25 


18.4± 1.8 


105.6±10.0 


134.3± 13.5 


5.1 


xlO" 


-9 


MPl 


644±28 


15.3± 3.5 


2.2±4.6 


6.4± 5.5 


1.9 


xlO" 


-5 


MP2 


586± 6 


6.2± 0.5 


-3.0±4.4 


-1.6± 4.7 


3.8 


xlO" 


-5 


MPS 


694± 8 


33.4± 2.5 


-1.7±4.5 


7.6± 7.0 


8.0 


xlO-^ 


-6 


MP4 


905±39 


23. 7± 4.8 


12.7±4.8 


20.9± 6.7 


9.3 


xlO" 


-8 


MP5 


739±15 


21. 5± 2.6 


0.5±4.5 


5.8± 5.9 


2.2 


xlO" 


-6 


MP6 


599± 7 


32.6± 2.6 


-1.1±4.5 


9.0± 6.8 


1.4 


xlO" 


-5 


MP14 


858±25 


15.0± 2.2 


11.2±4.8 


16.1± 5.8 


3.1 


xlO-^ 


-7 



Note. — Dust properties derived from the SED fitting of 11 WCs in six bands: 
J, H, F190N, Kg, 3.6/im and 4.5 fim. The second column presents the single dust 
temperature, T^ust (K)- The percentage of the total F190N flux attributable to dust 
is listed in the third column. The fourth column shows the uncorrected EW as listed 
in Table [21 while the fifth column presents the corrected value (after subtracting the 
dust contribution). The derived dust mass is listed in the sixth column. 



